Comment on "Comparison of the Cc and R3c space groups for the superlattice phase 

Of Pb(Zr .52Tio.48)0 3 " 
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The two-phase (space groups Cc and Cm) structural model recently proposed for the modeling 
of the neutron powder diffaction pattern collected on Pb(Zr a; Tii_ a ;)03 (PZT) sample with x ~ 0.52 
[R. Ranjan, A. K. Singh, Ragini, and D. Pandey. Phys. Rev. B 71, 092101 (2005)] is shown to 
neglect the hkl dependent (anisotropic) line broadening. The most serious problem with this model 
is that it assigns octahedral tilts to a wrong phase. Instead of correctly taking the anisotropic line 
broadening into account, this model used low symmetry phases to minimize the residuals during 
Rietveld refinement. The essential features of a model taking the hkl dependent line broadening into 
account are summarized. It has already been reported that once the anisotropic line broadening, 
revealed by high resolution neutron powder diffraction instrument, is correctly taken into account 
the model with R3c and Cm space group symmetries describes all Bragg peaks and their intensities 
well [J. Frantti, S. Eriksson, S. Hull, V. Lantto, H. Rundlof, and M. Kakihana. J. Phys.: Condens. 
Matter 15, 6031 (2003).]. It was further shown that the Cm + R3c model is consistent with the 
structural features observed at other compositions and temperatures, which is particularly important 
in the vicinity of the phase boundary. Also problems related to sample preparation and data 
collection are pointed out. 

PACS numbers: 77.84.Dy 61.12.Ld 61.50.Ah 81.30.Dz 



Recently, based on previously published neutron pow- 
der diffraction datai^, a model which used two mono- 
clinic phases (space group symmetries Cc and Cm) was 
proposed to describe the structure of Pb(Zr :c Tii_ ;I ;)03 
(PZT) with x » 0.52 3 . This model assigned a space 
group symmetries P4mm (at room temperature) or Cc 
(at low temperature) to the phase which was assigned to 
Cm symmetry in refs. and|J In addition, the phase, 
which traditionally has been assigned to the rhombohc- 
dral symmetry, was assigned to Cm symmetry. To un- 
derstand the problems related to such a model and the 
reason why the use of space group Cc was previously 
rejected 5 it is necessary to summarize the crucial role of 
the hkl dependent line broadening for Rietveld refine- 
ment in the case of these materials. There is a high risk 
that erroneous space group symmetry assignments result 
in once this line broadening is compensated for by reduc- 
ing the space group symmetry. For example, one can 
make the difference between measured and calculated 
diffraction intensities arbitrarily small by simply using 
space group PI and increasing the primitive cell size un- 
til the difference vanishes below the desired value. In our 
opinion, a structure model which reduces the symmetries 
to 'model' the hkl dependent line broadening takes the 
disorder into account in a wrong way, since in this way the 
disorder is assigned to be a periodical disturbance extend- 
ing through the crystal, obeying the space group symme- 
try. However, the /ifcZ-dependent line broadening is an 
inherent property of PZT powders, and is already seen 
in tetragonal Ti rich compositions. The same behaviour 
was observed in a closely related Pb(Hf 2 .Tii_ x )03 (PHT) 
system with 0.10 < x < 0.40, where the peak widths of 
the 001 reflections were twice as large as the widths of 



hOO reflections^. It must be emphasized that no signs of 
symmetry lowering from the P4mm symmetry was ob- 
served in this high- resolution neutron powder diffraction 
study. Light scattering experiments revealed that both 
in the case of PZT and PHT samples deviations from 
the average symmetry were observedi^. Although one 
could decrease the average symmetry to take these obser- 
vations into account, a more realistic model is to assume 
that the average symmetry over a length scale of a few 
hundred nanometers is P4mm, and deviations from this 
symmetry occur in a local scale (of a few unit cells) . The 
essential point here, too, is that these deviations are not 
periodical! It is also important to note that the Cm phase 
was identified by studying the peak split of the pseudo- 
cubic 110 reflections 4 . This peak split was beyond the 
resolution of the instrument used in ref. and there is 
no way to recover the lost information. It does not mat- 
ter how many 'new' refinements are carried out for the 
low-resolution data collected with insufficient counting 
time (note that neutron powder diffraction data in refs. 
[ijy and y are the same) . To solve the problem related 
to average symmetries it is necessary to carry out ex- 
periments with an appropriate instrument and carefully 
prepared sample. The role of the spatial composition 
variation (which in practice cannot be eliminated) in the 
vicinity of the phase boundary is an old and still a valid 
explanation for the two phase 'co-existences' and also 
partially explains the /ifcZ-dependent line broadening. It 
was also worth of trying to find a simple model which is 
consistent with other compositions and temperature s 5 ' 12 . 
These aspects are shortly reviewed below. 

a. Spatial composition variation. Traditionally, the 
two-phase 'co-existence', observed in room temperature 
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FIG. 1: The consequences of the spatial composition varia- 
tion (at fixed temperature and pressure) in the case of two 
samples with average compositions xi and X2 in the vicinity 
of (a) the phase boundary separating R3c and Cm phases (at 
X\) and (b) two phase boundaries separating R3c and Cc (at 
X2) and Cc and Cm phases (at -X3). 



PZT in the vicinity of the morphotropic phase boundary 
(MPB), has been explained by the composition variation. 
In the context of low temperature symmetries the role of 
spatial composition variation was discussed in refs. .5 
and 0. The crystal symmetry of Zr rich PZT ceram- 
ics is .R3(ii2ii2iii, except for the compositions with x m 1. 
Since spatial composition variation cannot be completely 
eliminated, there must exist two phases in the vicinity of 
MPB. To allow the existence of Cc phase necessitates 
that there should be a narrow region in the x — T plane 
were this phase is stable or metastable. Now, if we were 
to explain the existence of Cc phase, two phase bound- 
aries located somewhere between 0.52 < x < 0.54 should 
be assumed to exists, see Fig. Q](b). This in turn, once 
the spatial composition variation is taken into account, 
leads to three phase 'co-existence'. We preferred the sim- 
plest Cm+R3c model (corresponding to Fig. n( a ))> since 
it was able to explain the experimental observations in 
simplest terms. At low temperature (4 or 10 K£) 
the phase fraction of the Cm phase was monotonically 
decreasing with x increasing from 0.52 to 0.54, which 
in turn implies that two-phase 'co-existence' is predom- 
inantly due to the spatial composition variation. Other 
factors include stresses, which probably exist in the grain 
boundaries. It is worth to point out that the model pro- 
posed in ref. was ruled out also in a recent paper, see 
the last paragraph in ref. 0. After all, the only differ- 
ence between our model and the model proposed in ref. 
IT5I is that we used R3c symmetry (corresponding to two 
lattice parameters and four atomic coordinates), while 
they used Cc symmetry (corresponding to four lattice 
parameters and twelve atomic coordinates, if constraints 
are not used), a subgroup of R3c space group. Thus, 
the model proposed by us and the model proposed in 
ref. assigned octahedral tilts (evidenced by superlat- 
tice reflections) to the same phase. Octahedral tilts are 
crucial once the structural changes versus temperature, 
composition or pressure are studied. Somewhat surpris- 
ingly, the space group symmetries of the two phases were 
rather arbitrarily interchanged in refs. and IToL This 
model immediately runs into trouble, as the octahedral 
tilts are assigned to a wrong phase, see below. 



In addition, there have been reports claiming that di- 
electric measurements support Cc phase (for example, 
see ref. 0). In this context we note that the samples 
studied in refs. and ^3 contained a nonidentified im- 
purity phase(s), as was revealed by the peaks at around 
28 and 35 two-theta degrees (corresponding to d spacings 
3.18 A and 2.56 A, respectively). This implies that the 
composition was not well known and there is no way to 
avoid compositional and structural inhomogeneities. Al- 
though this non-perovskite phase was clearly observable, 
it was neglected in the Rietveld refinement model con- 
sidered in ref. [ToL which in turn results in an error in the 
structural parameters of the perovskite phase(s). This 
is related to the fact that the diffraction pattern shown 
in ref. [h] (reported to have x = 0.52) is reminiscent to 
our diffraction pattern with x — 0.53. In contrast, our 
diffraction pattern with x — 0.52 (see ref. ITlf) and the 
diffraction pattern of the x — 0.52 sample reported in 
refs. and were reminiscent. It should also be un- 
derstood that dielectric measurements do not suit for a 
space group determination. Although the dielectric mea- 
surements can provide valuable information in the case 
of single phase samples, particular care is necessary once 
multiphase samples are studied. Anomalies observed in 
a dielectric constant versus temperature measured from 
a two or three phase ceramic bulk samples (see ref. Il7|) 
do not provide reliable evidence that one of the phases 
has undergone a phase transition from Cm to Cc phase. 

b. Anisotropic line broadening. The anisotropic line 
broadening was ascribed to the spatial composition vari- 
ation resulting in 'microstrain' (i.e., Zr substitution for 
Ti creates local strains which in turn contributes to the 
hkl dependent line broadening)^. In addition, rather 
strong diffuse scattering is commonly observed between 
the hOO and OO/i Bragg reflections, which was assigned 
to Pb ions displaced toward (110) directions. This lat- 
ter feature is somewhat puzzling, as it adds intensity to 
certain diffraction peaks, which results in an asymmetric 
profile. For example, the larger and smaller d-spacings 
sides of the pseudo-cubic 200 and 002 reflections, respec- 
tively, gain intensity. In ref. Cm phase was practically 
used to model the intensity not only due to the crystalline 
phase, but also due to the diffuse scattering and due to 
the main crystalline phase. Despite the increased num- 
ber of refined parameters, the residuals were still rather 
high and the differences between different models com- 
pared were marginal (for instance, R exp was lowest for 
the R3c + Cm model which was rejected in ref. |3j). 

Although a reasonable structure refinement for 
medium resolution data collected on PZT samples with 
x < 0.50 could be obtained using a lineshape which ig- 
nores anisotropic line broadening, the situation was quite 
different for high resolution facilities, particularly once 
the data was collected on two phase samples (the present 
case). In such a case it was essential to use an appropri- 
ate profile function. Now, if the anisotropic line broaden- 
ing is neglected, the fit in the case of certain weak peaks 
becomes slightly worse (as was observed to be the case 
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FIG. 2: Left panel: Schematic illustration of the way how the 
neglect of anisotropic line broadening affects the position of 
the rhombohedral peaks: once the rightmost peak broadens, 
it pushes the rhombohedral peak (middle) towards smaller 
d-spacing. Righ panel: Observed and computed high res- 
olution neutron powder diffraction profile collected on PZT 
with x = 0.53 at 4 K. Pseudo-cubic 200 reflection region is 
shown. Note that the 002m peak (Cm phase) at around 2.06 
A is significantly broader than the 220m (Cm phase) peak at 
around 2.02 A. Now, the position of the rhombohedral 204_r 
peak depends on the way the line broadening is taken into 
account. 

of weak superlattice reflections, which are less weighted 
in the refinements). The reason for this is illustrated 
in Fig. [21 in order to improve the fit corresponding to 
the strong pseudo-cubic 200 reflections, during the re- 
finement the position of R3c reflection is shifted toward 
higher d-spacings, which in turn resulted in a small shift 
of the weak peaks at 2.44 and 1.06 A (in a case of the 
data shown in ref. 3j the latter peak was almost at the 
level of noise and the Cc + Cm model assigns more re- 
flections than there are data points in this region). The 
anisotropic line broadening was not limited to the d- 
spacing at around 2 A, and similar mechanism was seen 
at other d-spacings. 

c. Previous high resolution neutron powder diffrac- 
tion results. The anisotropic line broadening was the 
reason to carry out a subsequent study using a high reso- 
lution neutron powder diffraction instrument (which re- 
veals the anisotropic line broadening particularly clearly 
and allows a distinction between average symmetry low- 
ering and a defect related anisotropic line broadening 
to be done, in contrast to the low- or medium resolu- 
tion constant wavelength neutron diffraction facilities) 12 . 



In addition, we increased Zr content, x — 0.54^ (al- 
though it definitely was possible to fit all reflections us- 
ing the Cm + R3c model also in the case of x — 0.52 and 
x = 0.53 samples). This allowed a more reliable refine- 
ment and symmetry identification to be done by studying 
the changes in phase fractions versus temperature. Im- 
portantly, for these compositions it has been found that 
oxygen octahedra tilts increase with increasing x (ref. 
UM and decreasing temperature (this feature is discussed 
in refs. and@). We also note that previously Cc was 
proposed to be space group symmetry corresponding to 
high isotropic pressure^, whereas our work was concen- 
trating on the determination of space group symmetries 
versus composition and temperature at ambient pressure. 
Now, PZT sample with x = 0.54 provided a test for clar- 
ifying which phase is the preferred one at low tempera- 
ture. To model the peak profiles, GSASi^ lineshape 4 by 
Stephens^ was used. As was seen in ref. 0, once the 
crystal contracts with decreasing temperature the i?3c 
phase was favoured (and Cm phase did not transform to 
Cc phase, as one could expect from the model proposed 
in ref. 3j. Thus, when the oxygen octahedra had almost 
no chance to contract, they were tilt towards opposite 
directions along the pseudo-cubic 111 axis (tilt system 
a~a~a~). This tilt results in a decrease of the volume 
of cuboctahedra around Pb ions. Changes evidencing 
Cc symmetry (and the corresponding tilt system) were 
not observed. Indeed, the phase fraction of Cm phase 
(which was assigned to Cc symmetry in ref. \Q) decreased 
with decreasing temperature, whereas i?3c phase fraction 
significantly increased with decreasing temperature. This 
can be confirmed even by a naked eye, although the re- 
sults of refinements were also given. The intensity of 
the peak at around 1.06 A was increasing with decreas- 
ing temperature, which further confirmed that its origin is 
the phase assigned to R3c phase. Also the peak at around 
1.06 A was well fit by R3c symmetry. 

As a summary, we do not find support for Cm + Cc 
model proposed in ref. 0. Instead, we still consider that 
Cm+R3c provides a sufficient low temperature structural 
model for PZT with composition in the vicinity of the 
morphotropic phase boundary. 
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